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ABSTRACT 

The  model  of  a  Surface  Effect  Ship  was  refined  to  include 
simplified  propulsion  dynamics,  negative  drag  characteris- 
tics, sea  state  effects  and  an  autopilot  for  speed  control. 
These  design  modifications  were  introduced  into  a  real  time, 
man  controlled  simulation  of  a  3000  ton  Surface  Effect  Ship 
(3K-SES)  in  5  degrees  of  freedom  (RTS5D)  and  results  were 
compared  with  a  Data  Base  Program  (DBSIM5D)  based  on  towing 
tank  data  scaled  up  to  model  the  ship. 

Hardware  and  Software  design  changes  were  incorporated 
into  the  RTS5D  model  to  provide  a  more  accurate  approxima- 
tion of  real  time,  a  faster  computer  iteration  time,  and  a 
broach  condition  warning  if  the  operator  exceeded  certain 
thrust  vectoring  limits. 
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NOMENCLATURE 


A      Average  sidewall  wetted  area,        ft2 


ws 


starboard  side 


9 

A      Average  sidewall  wetted  area,  •       ft" 


wp 


port  side 


A31    Added  mass  coefficient  in  roll       ft-slug 
force  equation 

A33    Added  mass  coefficient  in  pitch      ft-slug 
force  equation 

2 
A,    Average  wetted  sidewall  area         ft 


w2 


of  the  bow 


2 
A  ,     Average  wetted  sidewall  area         ft 


wl 


of  the  stern 


A22    Added  mass  coefficient  in  yaw        slug  s   lbf/ 
force  equation  ft 

8      Sideslip  angle  rad 

2      2 
C^y    Coefficients  drag   in  x-direction    lb^s  lb  /ft"" 

DA(K)  t  m 

2      2 
C^v    Coefficient  drag  in  y-direction      lb,.s**lb  /ft 
DY  °  t    m 

Cn7_   Sidewall  roll  moment  lumped 
parameter  coefficient 

Cnp    Bow  pitch  force  lumped  parameter 
coefficient 

F      Sidewall  roll  force 
sw 

F      Sidewall  starboard  buoyancy  force 
ss  J 

F  Sidewall  port  buoyancy  force 

F  Stern  buoyancy  force 

F?  Bow  seal  pitching  force 

F  Bow  buoyancy  force 

g  Gravitational  acceleration 


2    2 
lb£sVft^ 

non- 

dimensiona 

1 

lbs 

lbs 

lbs 

lbs 

lbs 

lbs 

ft/s2 

2 
l„      Moment  of  inertia  about  x-axis      slug  ft 

2 
Iy      Moment  of  intertia  about  y-axis     slug  ft 

2 
I7      Moment  of  inertia  about  z-axis      slug  ft 

K       Summation  of  moments  about  x-axis    lbs-ft 

I  Length  of  sidewall  ft 

I,  Actual  draft  of  port  sidewall       ft 

I  x-direction  displacement  of  hull    ft 

drag  centroid 

lA  Average  draft  of  SES  sidewall       ft 

dl 
i.      Average  draft  of  bow  seal  ft 

I  Length  from  center  of  gravity  to 

stern 

l  Average  draft  of  widewall  ft 

xi 

£,,      Pitch  moment  lever  arm  for  bow      ft 
sidewall  buoyance  force 

l-  Pitch  moment  lever  arm  for  bow      ft 

seal  force 

m  Mass  of  the  rigid  ship 

M  Summation  of  moments  about  y-axis 

N  Summation  of  moments  about  z-axis 

P,  Plenum  pressure 

D 

p  Lumped  drag  centroid  point 
p*  Lumped  drag  centroid  point 
p       Roll  acceleration 


ft 

lbs- 

ft 

lbs- 

ft 

lbs- 

ft2 

non- 

dimensiona 

1 

non- 

dimensiona 

1 

rad/ 

2 
sec1" 
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q      Pitch  rate  rad/sec 

q      Pitch  acceleration  rad/sec2 

r      Yaw  rate  rad/sec 

r      Yaw  acceleration  rad/sec 

s^     Turning  moment  lever  arm  of  no.  1    ft 
engine 

s2     Turning  moment  lever  arm  of  no.  2    ft 
engine 

s-     Turning  moment  lever  arm  of  no.  3    ft 
engine 

s.     Turning  moment  lever  arm  of  no.  4    ft 
engine 

T_     Total  thrust  magnitude  on  no .  1      lbs 
engine 

Tg     Total  thrust  magnitude  on  no .  2      lbs 


-©  ■ 


engine 


T,     Total  thrust  magnitude  on  no .  3      lbs 
engine 

T-Q    Total  thrust  magnitude  on  no .  4      lbs 
engine 

Tr     Total  forward  thrust  vector  of       lbs 

fOrW      rr 

effectors 

T  .  ,   Total  side  thrust  vector  of  lbs 

side    Cr 

effectors 

T      Total  turning  moment  generated  by    lbs  ft 
^     effectors 

u      Velocity  in  x-direction  (surge)      ft/sec 

2 
u      Acceleration  in  x-direction  ft/sec 

v      Velocity  in  y-direction  (sway)       ft/sec 

2 
v      Acceleration  in  y-direction  ft/sec 

v      Total  velocity  ft/sec 
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Width  of  bow  seal  ft 

X      Summation  of  forces  in  x-direction   lbs 


w 
e 


Xo  XNAV  coordinate  of  SES  ft 

X  XATA„  velocity  of  SES  ft/sec 

o  NAV 

Y  Summation  of  forces  in  y-direction   lbs 

Y  YXTA,r  coordinates  of  SES  ft 

o  NAV 

Y  YXT.,,  velocity  of  SES  ft/sec 

o  NAV        J 

6  Effector  angle  commanded  rad 

5y  Effector  angle  of  no.  1  nozzle       rad 

<5„  Effector  angle  of  no.  2  nozzle       rad 

6q  Effector  angle  of  no.  3  nozzle       rad 

5  Effector  angle  of  no.  4  nozzle       rad 

i>  Heading  angle  of  SES  rad 

4>  Roll  angle  of  SES  rad 

9  Pitch  angle  of  SES  rad 

p  Density  of  sea  water  lbm/ft' 
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I.   INTRODUCTION 

Rising  fuel  costs,  critical  manpower  shortages  and  an 
increasing  need  for  rapid  deployment  of  forces  in  defense 
of  United  States'  interests  throughout  the  world  mandate 
the  development  of  fuel  efficient,  high  speed,  minimum  man- 
ned  ships  for  the  Navy.   The  Captured  Air  Buddie  (CAB)  Sur- 
face Effect  Ship  (SES)  holds  great  promise  in  providing  an 
answer  to  these  problems. 

A  ship  with  minimum  manning  traveling  at  the  high  speeds 
envisioned  for  the  SES  would  require  highly  skilled  operator 
personnel.   It  would  be  necessary  to  provide  these  personnel 
with  real  time  operational  training  similar  to  that  which 
pilots  receive  in  aircraft  simulators.   This  would  test  the 
man  in  responding  to  failure  situations  that  may  be  encoun- 
tered during  actual  operation  and  train  personnel  in  control 
techniques  that  are  unique  to  SES  type  vehicles.   A  Real  Time 
Simulator  could  also  be  used  as  a  means  to  test  hardware  used 
in  SES  craft  and  provide  a  design  tool  for  the  development 
of  future  modifications. 

A  simplified,  non-linear  five  degree  of  freedom  Real  Time 
Simulator  (RTS5D)  for  the  3K-SES  has  been  developed  by  T.  S. 
Nelson  [Ref.  1]  such  that  a  continuously  observable  real  time 
solution  is  interfaced  with  man-generated  control.   Certain 
"man-in-the-loop"  experiments  have  been  conducted  using  the 
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RTS5D  that  have  demonstrated  limitations  in  that  model.   It 
was  necessary  to:  1)  refine  ship  dynamics;  2)  introduce  sim- 
plified propulsion  dynamics;  3)  constrain  thrust  angle  in 
order  to  remain  within  the  bounds  of  programmed  dynamic  res- 
ponse (avoid  broaching) ;  4)  refine  the  software  to  provide 
faster  computer  iteration  time;  and  5)  change  the  hardware 
to  provide  the  operator  with  speed  control  as  well  as  turning 
control . 

Accordingly,  changes  have  been  made  in  the  RTS5D  model 
to  more  accurately  reflect  the  nonlinear  forward  drag 
characteristics  of  the  SES.   A  speed  controller  has  been 
designed  to  allow  operation  of  the  3K-SES  within  a  speed 
range  of  40-60  knots  and  limits  have  been  placed  on  the 
thruster  deflection  angle  in  order  to  prevent  broaching,  a 
condition  which  is  not  allowed  under  the  existing  dynamic 
model.   Software  refinements  have  been  incorporated  to  im- 
prove computer  iteration  time. 
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II.   EQUATIONS  OF  MOTION 

A.   COORDINATE  SYSTEMS  AND  ASSUMPTIONS 

The  RTS5D  model  of  the  3K-SES  used  simplified  equations  of 
motion.   For  completeness  of  comparison  the  development  of 
these  equations  is  repeated  here.   The  rational  for  requiring 
simplification  was  to  start  with  differential  equations  of 
motion  that  could  be  utilized  effectively  in  a  real  time  simu- 
lation.  To  satisfy  this  requirement  a  model  was  selected  that 
would  use  point  source  of  drag  forces  consisting  of  only  a 
few  lumped-coefficient  terms.   The  equations  used  are  based 
upon  those  developed  in  Ref.  2  which  follow  this  concept. 
Additional  assumptions  used  in  the  simplified  5  DOF  equations 
were  as  follows: 

1)  All  accelerations  are  measured  at  the  center  of 
gravity. 

2)  Vehicle  is  "free- to-heave" . 

3)  All  cross  coupled  moments  of  inertia  are  zero. 

4)  All  cross  products  of  angular  velocities  in  force 
equations  are  zero. 

5)  All  roll  and  pitch  angles  used  in  force  equations 
are  subject  to  small  angle  approximations. 

6)  Calm  water  conditions. 

7)  Mass  and  mass  distribution  of  vehicle  are  constant. 

8)  Effect  of  changes  of  aerodynamic  force  are  neglected. 
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Utilizing  the  above  assumptions,  the  simplified  5  DOF 
equations  of  motion  are  developed  using  the  coordinate  sys 
terns  shown  in  Fig.  1  and  Fig.  2,  the  force  vector  diagrams 
shown  in  Fig.  3,  Fig.  4A,  Fig.  4B,  and  the  following 
equations  from  Newton's  laws  of  motion: 


SURGE 

SWAY 

YAW 

PITCH 

ROLL 


m(u  -  vr) 
m(v  -  ur) 


h  r 


Ix  P 


=  X 

=  Y 

=  N 

=  M 

=  K 


where 
m 
u 
v 
r 

P 

q 

i. 

i. 


x 


X 
Y 
N 

M 
K 


mass  of  the  rigid  ship 

velocity  in  x-direction  (SURGE) 

velocity  in  y-direction  (SWAY) 

angular  velocity  about  the  z-axis 

angular  velocity  about  the  x-axis 

angular  velocity  about  the  y-axis 

moment  of  inertia  about  the  z-axis 

moment  of  inertia  about  the  y-axis 

moment  of  inertia  about  the  x-axis 

summation  of  forces  in  x-direction 

summation  of  forces  in  y-direction 

summation  of  moments  about  the 
z -  ax  i  s 

summation  of  moments  about  the  y-axis 

summation  of  moments  about  the  x-axis 


'  slugs 

'  ft/sec 

*  ft/sec 

'  deg/sec 

'  deg/sec 

'  deg/sec 

2 


•  lb  -ft 

m 

•  lb  -f 

m 

»  lb  -ft' 
m 


2 


lb 


lb 


•  ft-lb 


'  ft-lb 


'  ft-lb 
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Additionally  the  following  navigation  relationships  are 
utilized  (see  Fig.  1). 


X    =  u  cosijj   -  v  sin^ 
o 


Y  =  u   sinili      +   v  cosii; 

o 


2  2  2 

V  =     U        +     V 

s 


6  =    tan    ~1    (-v/u) 


where : 

ty   =  heading  angle 

6  =  drift  angle 

5  =  thrust  vector  angle 

Also, 


X 

0 

= 

/X  dt 
0 

Y 

0 

= 

/Y      dt 

0 

r 

= 

It   dt 

P 

= 

/p    dt 

q 

= 

/q    dt 

u  =  /u  dt 

v  =  /v  dt 

ty    =  /r  dt 

<f>  =  /p  dt 

9  =  /q  dt 
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Note 


1) 
2) 


4AV 


$-0   when  vehicle  heads  paralled  to  +  X 


axis . 


NAV 


Vehicle  in  figure  above  shown  with  negative 
thrust  vector  angle  6,  positive  \\i ,    negative 
8,  negative  sway  velocity  v,  i.e.,  in  a  right 
turn 


Figure  1 
Definition  Of  Coordinate  System  (Part  I) 
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Note:   Direction  for  positive  velocity  and 
angular  rates  are  shown 


Figure  2 
Definition  Of  Coordinate  System  C?art  II) 
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c 


+  r 


BOW  ! 


^CDYv|v 


PORT 
SIDEWALL 


^CDXU 


*CDXU 


w 


>5CDYV  I 


rpl 


STARBOARD 
SIDEWALL 


rTiosin<5io 


T10COS510 


xNote:   1)  p  and  p'  are  equivalent  point  force  centroids 

2)  5_,  5g,  5  ,  5  _  are  thrust  vector  angles  and 

are  not  required  to  be  equal. 

3)  T_,  T„,  Tg ,  T. 0  are  thrust  magnitudes  and  are 

not  required  to  be  equal. 

4)  Force  directions  shown  are  for  a  right  turn 

Figure  3 
Surface  Effect  Ship  CTop  View) 
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Note:   1)  Right  Turn  Forces  Acting,  i.e.,  Sj,    5g,  5g, 
5, n  are  negative 

2>  TSIDE  =  T7sinV+T8sin68  +  T95in59  +  T103inS10 


Figure  4A 
Surface  Effect  Ship  CStern  View) 
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+  z 


\ 


+  a 


FORW 


Note:   TFORW=T7COS57+T8COS58+T9COS69+T10COS510 
(see  Fig.  3) 


Figure  4B 
Surface  Effect  Ship  CSide  View) 
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B.   FORCES  AND  MOMENTS 

1.  Surge  Forces  (see  Fig.  4B) 

The  forward  acceleration  was  determined  to  be  genera- 
ted by  the  summation  of  thrust  vectored  in  the  forward  direc- 
tion  Tforw'   The  resultant  acceleration  of  the  SES  was 
assumed  to  be  opposed  by  a  retarding  force  exhibiting  a 
"velocity-squared"  characteristic.   Additionally,  Newton's 
laws  of  motion  specified  a  "v.r"  product  that  contributed 
to  the  drag  component  of  the  equation  when  the  vehicle  was  in 
a  turn  (contrifugal  force  reaction  term).   The  simplified 
equation  of  motion  for  the  surge  acceleration  follows  from 
a  summation  of  forces  in  the  u  direction  as  depicted  in 
Fig.  4B. 

u  =  (Tforw/m)  "  (CDXu2/m)  +  vr 

2.  Sway  Forces  (see  Fig.  4A) 

The  side  acceleration  of  the  vehicle  was  analyzed 
to  be  the  result  of  the  summation  of  the  thrust  vectored 
parallel  to  the  vehicle's  y-axis  and  a  retardation  "velocity 
squared"  phenomena  such  as  that  produced  by  the  cross  flow 
drag  term  of  the  sidewall  as  described  in  Ref.  3.   This 
retardation  force  is  augmented  by  a  "u-r"  product  as  speci- 
fied by  Newton's  Law. 


v  =  Tside/m  "  CDY  vlv'/m  "  ur 
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2 
Note  the  requirement  to  force  the  Cnyv  /m  term  to 

maintain  the  sign  of  the  sway  velocity  such  that  both  left 

and  right  turns  may  be  computed. 

3.   Yaw  Moments  (see  Fig.  3) 

The  yaw  acceleration  was  a  turning  moment  summation 

which  was  defined  to  be  influenced  by  three  primary  forces 

operating  over  moment  arms.   It  is  of  interest  to  note  that 

the  third  term  is  an  "added  mass"  term  which  is  considered 

to  operate  on  the  same  pressure  point  as  the  C  yv|v|  term. 

The  term  T    is  a  summation  of  moments  generated  by  the 
yaw  *  J 

four  thrusters  as  shown  in  Fig.  3. 


T    =  s.T7  cos5-  +  sTT0  cos  50  -  s..Tn  cos  5n 
yaw    47     7     38      8     19      9 


"  S2T10  C0S  610  •  *oT7  Sin  67  "  *oT8  Sin  68 


*oT9  sin  6g  -  £QT10  sin  6±Q 


r  =  CDY  zw   v|v|/Iz  +  A22uv*w/Iz  +  TyAW/I2 

4.   Pitch  Moments  (see  Fig.  4B) 

The  pitch  acceleration  was  assumed  to  be  the  summa- 
tion of  moments  generated  by  forward  thrust,  Tforw>  the 
buoyancy  of  the  sidewalls  of  the  SES,  F^  and  F-   and  a 
vertical  force  generated  at  the  bow  of  the  vehicle,  F_. 
This  vertical  bow  force  was  defined  to  be  a  lumped  parameter 
term  which  modeled  the  reaction  force  due  to  plenum  pressure 
acting  against  the  bow  seal. 
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l^        =  average  draft  of  bow  seal 
A    =  r  £   +  ((£  tane)/2)£ 

w,      XX,     ^  x      J       J    X 


A  =  £  £    -  (O   tane)/2£ 
^2            xx,      x         x 

F,  =  A    l  pg 

1      W]_   x2  pS 


F~    =  A    £    pg 

2      ™2   x2 


F3    =  CDpPbWe(£d  "  '31  tan  9) 


£  ,  =  average  draft  of  sidewall 

p,  =  plenum  pressure 

A  =  average  wetted  sidewall  area  of  the  stern 

wi 

£„,,  =  width  of  one  sidewall 
*2 


w  =  width  of  bow  seal 
e 

A  =  average  wetted  sidewall  area  of  the  bow 

w2 

A-.,  =  added  mass  coefficient 
3j> 

p  =  water  density 

g  =  32.3  ft/s2 

£..,  =  lever  arm  of  bow  seal 

lT  =  lever  arm  of  buoyancy  force 


q  =  (Tr   £„  +  F7^-i  +  F^  -  cnv  u2  l, 

n    v  forw   p     3  jl     23     Da      z 


"  Fl*3  "  A33Uq)/IY 
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Note  the  added  mass  term  A-^uq  which  was  required 
to  provide  damping  to  the  pitch  moment  and  is  specified 
in  Ref.  2. 

It  was  found  that  the  response  of  the  SES  to  a 
pitch  perturbation  without  the  added  mass  component  was 
undamped  and  approximately  sinusoidal.   Unlike  the  yaw 
equation  where  the  added  mass  term  was  a  small  contributor 
to  damping,  the  pitch  added  mass  term  was  found  to  be  of 
significant  importance  in  modeling  the  known  pitch  motion. 
5.   Roll  Moments  (see  Fig.  4A) 

The  roll  acceleration  equation  was  analyzed  to  be  a 

summation  of  moments  generated  by  buoyancy  forces  of  the 

port  and  starboard  sidewalls,  F    and  F   ,  the  thrust 
r  sp      ss 

vectored  parallel  to  the  y-axis  of  the  vehicle,  T  .,  ,  a  side 
r  J  side 

force  C  y  v|v|,  and  a  lumped  coefficient  vertical  force, 

F   .   The  vertical  force  F   was  defined  as  the  force 

sw  sw 

generated  in  a  turn  due  to  the  sidewall  curvature  (dead 

rise  angle)  acting  against  the  cross  flow  of  water. 

F     =  pg  A   £dp 
sp      6   wp   r 

F     =  pg  A   £ds 
ss      &   ws 

F     =  Cn7D  v|v| 
sw      DZP 

i,    =  average  draft  of  SES  sidewall 

H,n   =  lA  *  r  tan  4> 

dp      a-.  -   sw 

l.  =  I  *      +   I        tan  <$> 

ds      d.     sw 
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A    =  average  wetted  area  port 

A    =  average  wetted  area  starboard 


p=((F    -F  )l        -  T  .,  I      +  C„„  vlvl  % 
v        v  ^  sp    ss;  sw    side  p    DY   '  '   z 


-  F  I        -   A,n  up)/I__ 
sw  sw    31  VJI    X 

Again  note  the  added  mass  term  A.,_up/I„  which  was 
found  to  be  essential  in  modeling  the  damping  phenomena  of 
the  roll  motion. 

In  summary,  the  simplified  5  DOF  equations  of  motion 
used  in  RTS5D  for  the  3K  TON  SES  are: 

SURGE   u  =  Tforw/m  -  CDX  u2/m  +  vr  (Fig.  4B) 

SWAY    v  =  T  .,  /m  -  CL.V  v|v|/m  =  ur        (Fig.  4A) 

s lde      ui 

XM     r  =  Tyaw/Iz  +  CDY  v|v|  £w/Iz  *  A22uv.w/Iz 

(Fig.  3) 

PITCH    q  =  (T,    £   +  F-^-zi  +  F7*t  "  Cr,Y  u~Zz 

n      forw  p     :>  31     L    3     Ua 

-  ?1lz    -  A34uq)/Iy  (Fig.  4B) 

mk  P  -  CCFsp  -  Fss).sw  -  Tsideip    ♦  CDy  v|vU, 

"  Fsw*sw  "  A33  UP)/XX      fFiS-  4A) 
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C.   PARAMETER  IDENTIFICATION 

These  equations  are  an  extension  of  the  3  DOF  flat  turn 

SES  model  developed  by  Gerba  and  Thaler  in  Ref.  4.   The 

identification  of  craft  parameters  <!_.„,  C_..,  and  I      is 

r  DX    DY       w 

described  in  Ref.  4  and  repeated  here  for  completeness. 

The  surge  drag  coefficient  C^  is  determined  by  selecting 

a  steady  state  turn  condition,  where  T~        ,    u,  v,  r,  and  m 
7  '        forw'   '      * 

are  known. 

2 

n   Tforw   CDXU 

u  =  0  =  -  +  vr 

m       m 

from  which 

T  j.          +  mvr 
r    _   forw 

LDX        2 
u 

The  sway  drag  coefficient  CQY  is  determined  by  the  same 
steady  state  turn  condition  which  requires 

;  =  o  =  -^^  -  ^x  v|v|  .  ur 

m      m    '  ' 

where 

T.j   -  mur 
r        =   side 
DY       vjT| 

The  sway  drag  moment  arm  follows  utilizing  the  yaw  accel- 
eration equation 

•    n   Vw  ,  CDYY'vl'w  .  A2  2  UVV 

r  =  °  =    i7"    — r7 —    — U 
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which  yields 


yaw 
w  - 


CDY  v|v|  ♦  A22uv 


The  addition  of  the  pitch  and  roll  equations  introduces 

additional  lumped  coefficients  CDp  and  CDZp.   These  are 

easily  solved  using  known  constants  l    ,  l    ,  C^,  I    ,  £ 

w'   p'   DY'   z*   sw' 

w^>  Pu>  1a    and  steady  state  values  of  F   ,  F   ,  F,  ,  Fn, 
e,rD'd  J  sp   ss '   1    2 ' 

T^^^,,  T  .,  ,  u,  v,  and  r.   It  is  significant  to  note  that 
rorw    side  ° 

in  a  steady  state  condition  the  angular  roll  rate,  p,  and 
angular  pitch  rate,  q,  are  both  required  to  equal  zero. 
This  is  in  contrast  to  the  yaw  rate,  r,  where  a  finite 
steady  state  value  is  desired- in  a  turn.   Thus  the  added 
mass  terms  of  A,,uq  and  A_,  up  in  the  pitch  and  roll  equa- 
tions are  not  utilized  in  the  determination  of  craft  para- 
meters; their  function  is  strictly  confined  to  damping  of 
their  respective  accelerations.   Therefore,  utilizing  the 
pitch  equation  with  known  steady  state  turn  values. 

«  "  °  "  (Tfon/p  +  CDpPbwe(*d  "  l31  tan  6)131 

+  Vs  -  CDXu2*2  "  Vs^Y 
yields 

-  CDXAz    1V3    -    Vs    -    Tforw    S 

DP  PbV4d    "    *31    tan    8Jl31 
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The  roll  equation  under  steady  state  turn  conditions 
yields 


P  ■  0  -  C(Fsp  -  Fss)*w  -  Tside  ,p  *  CDY  v|v|  *z 


CDZP  VM  ^w^z 


from  which 


Cn7D  =((F    -F   )A   T-  .,  £   +  C^v  v   |v  \%    )/ 
DZP      sp    ss^  w  side  p    DY   ss1  ss1  z 


where 


^  ss '  ss '  sw; 


F  =  port  sidewall  buoyancy  in  pounds      (Fig.  4A) 

F  (  =  starboard  sidewall  buoyancy  in  pounds  (Fig.  4A) 

Ts . ,  =  thrust  vectored  parallel  to  vehicles  y-axis 

in  pounds  (Fig.  3) 

v  =  sway  steady  state  velocity  in  feet/sec 

C_Y  =  sway  drag  coefficient  for  steady  state  condition 

l  =  sway  drag  moment  in  feet              (Fig.  3) 

l  =  effector  thrust  moment  arm  in  feet     (Fig.  4A) 

l  =  surge  drag  moment  arm  in  feet         (Fig.  4A) 

l  =  sidewall  buoyancy  force  moment  arm 

sw  in  feet                             (Fig.  4A) 
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III.   RTS5D  VALIDATION 

A.  DBSIM5D  BENCHMARK 

Validation  of  the  RTS5D  was  attempted  using  the  data 
base  program  (DBSIM5D)  [Ref.  5]  as  a  benchmark.   However, 
new  guidance  in  the  proper  use  of  the  data  base  program  con- 
tained in  Ref.  6  revealed  the  following  errors  in  the  use  of 
the  BBSIM5D.   1)   Initial  conditions  had  been  improperly  set 
and  2)   Data  used  for  validation  for  the  56  knot  run  at 
effector  deflection  angles  greater  than  10°  was  not  correct 
because  of  the  thrust  inlet  broach  conditions  which  existed 
using  the  DBSIM5D.   These  conditions  were  not  accounted  for 
in  the  RTS5D  model.   The  noted  discrepancies  in  the  opera- 
tion of  the  DBSIM5D  necessitated  a  revalidation  of  the 
RTS5D. 

B.  VALIDATION  RESULTS 

Using  the  information  contained  in  Ref.  6,  the  DBSIM5D 
program  was  again  used  as  a  benchmark  to  validate  the  RTS5D 
simulation.   A  DBSIM5D  sequence  of  runs  with  initial  forward 
velocities  equal  to  40  knots,  50  knots  and  60  knots  with 
step  effector  angles  of  5°  and  10°  (at  40  knots),  5°,  10° 
and  15°  (at  50  knots),  and  5°  (at  60  knots)  were  compared 
to  identical  runs  of  the  RTS5D.   Thrust  effector  angles 
were  kept  below  a  maximum  that  would  have  caused  a  broach 
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condition.   Broach  conditions  are  explained  in  Chapter  IV. 
The  first  peak  overshoot  value  and  "Quasi  Steady  State" 
values  for  the  variables  u  (surge) ,  v  (sway) ,  r  (yaw  rate) , 
4>  (roll),  and  e  (pitch)  were  used  as  a  basis  for  comparison. 
The  results  are  shown  in  Tables  I-III.   Significant  error 
in  forward  velocity  occurred  in  the  50  knot  test.   Also,  it 
was  noted  that  in  40  knot  test  with  5°  thruster  deflection 
angle  that  the  DBSIM5D  actually  reached  a  "Quasi  Steady 
State"  forward  velocity  greater  than  the  initial  value  upon 
entering  the  turn.   These  conditions  led  to  the  development 
of  the  model  of  the  RTS5D  with  negative  drag  characteristics 
in  sea  state.   Differences  in  the  other  measured  parameters 
were  addressed  by  parameter  adjustment  on  the  new  model  in 
order  to  more  closely  match  values  obtained  with  the  DBSIM5D 
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TABLE  I 

RTS5D  and  DBSIM5D  Performance  Test  at  40  Knot: 

effector 
angle  DBSIM5D  RTS5D 


l=>  ' 


<j>                      in  degrees  <|> 

1st  pk  Qss                       1st  pk  Qss 

5°                 1.01  .742                      .17  .12 

10°                   1.46  1.39                        .31  .24 


$         in  degrees  <j> 

Qss  1st  pk  Qss 

1.40  1.27  1.21 

1.48  1.28  1.20 

u         in  ft/sec  u 

Qss  1st  pk  Qss 

70.9  n/a  67.03 

66.3  n/a  65.98 

v         in  ft/sec  v 

Qss  1st  pk  Qss 

1.37  2.84  2.42 

3.26  3.83  3.41 

r         in  degrees/sec  r 

Qss  1st  pk  Qss 

.403  .58  .33 

.959  .98  .66 


1st  pk 

5° 

1.44 

10° 

1.48 

1st  pk 

5° 

67.4 

10° 

n/a 

1st  pk 

5° 

1.42 

10° 

3.20 

1st  pk 

5° 

.396 

10° 

.953 
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TABLE  II 

RTS5D  and  DBSIM5D  Performance  Test  at  50  Knot: 


effector 
angle 

DBSIM5D 

RTS5D 

1st  pk 

4> 

in  degrees 
Qss 

1st  pk 

Qss 

5° 

.424 

.446 

.27 

.19 

10° 

.757 

.886 

.49 

.38 

15° 

1.02 

1.4 

.72 

.57 

1st  pk 

9 

in  deg' 
Qss 

rees 

9 
1st  pk 

Qss 

5° 

.862 

.892 

1.16 

1.15 

10° 

.943 

1.13 

1.15 

1.16 

15° 

1.12 

1.35 

1.22 

1.17 

1st  pk 

u 

in  ft/ 
Qss 

sec 

u 
1st  pk 

Qss 

5° 

n/a 

83.6 

n/a 

85.8 

10° 

n/a 

78.7 

n/a 

82.51 

15° 

n/a 

69.2 

n/a 

80.58 

1st  pk 

V 

in  ft/ 
Qss 

sec 

V 

1st  pk 

Qss 

5° 

1.31 

1.30 

3.55 

3.02 

10° 

2.72 

2.70 

4.79 

4.26 

15° 

4.05 

4.33 

5.69 

5.207 

1st  pk 

r 

in  deg 
Qss 

rees/ 

sec    r 

1st  pk 

Qss 

5° 

.368 

.365 

.73 

.41 

10° 

.864 

.870 

1.22 

.82 

15° 

1.35 

1.44 

1.66 

1.26 

34 


TABLE  III 

RTS5D  and 

DBSIM5D  Performance 

Test  at  60  Knots 

effector 

angle 

DBSIM5D 

RTS5D 

* 

in  degrees 

4> 

3° 

1st  pk 
.141 

9 

Qss 
.200 

in  degrees 

1st  pk 
.39 

9 

Qss 

.274 

5° 

1st  pk 

.521 

U 

Qss 
.583 

in  ft/sec 

1st  pk 
1.14 

u 

Qss 
1.07 

5° 

1st  pk 
n/a 

V 

Qss 
97.8 

in  ft/sec 

1st  pk 
n/a 

V 

Qss 
100.29 

5° 

1st  pk 
2.03 

Qss 
1.84 

1st  pk 
4.26 

Qss 
3.62 

r 

in  degrees/sec 

r 

5° 

1st  pk 
.552 

Qss 
.494 

1st  pk 
.87 

Qss 
.487 
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IV.   MODEL  DEVELOPMENT 

A.   DRAG  CHARACTERISTICS 

1 .  Surge  Equation  Model 

The  RTS5D  model  assumed  a  velocity  squared  drag  re- 
tarding characteristic  for  the  3K-SES  which  led  to  the  equa- 
tion for  surge  acceleration  that  was  presented  in  Chapter  II. 
The  drag  characteristics  from  Ref.  7  for  the  3K-SES  that  is 
simulated  by  the  DBSIM5D  are  shown  in  Fig.  5.   The  data  base 
program  was  run  using  parameters  that  caused  it  to  exhibit 
a  retarding  force  corresponding  to  the  curve  at  sea  state 
four,  therefore  this  curve  was  used  to  remodel  the  RTS5D. 
The  new  simplified  equation  of  motion  for  surge  acceleration 
became 

u  =  (Tforw/ml  "  (FSD/m)  +  v  •  r 

where  FSD  is  full  scale  drag  which  is  defined  by  the  curve 

in  Fig.  5.   T^-.,  is  the  sumation  of  the  four  thrusts  assumed 
rOKW 

in  this  model  to  be  applied  at  the  same  point  along  the 
center  line  of  the  craft. 

2.  Implementation 

FSD  was  implemented  into  the  RTS5D  model  by  dividing 
the  curve  into  four  separate  regions  and  defining  FSD  within 
each  region.   In  region  I  the  retarding  force  was  approxi- 
mated by  a  normalized  velocity  cubed  term.   In  regions  II 
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and  III  straight  line  approximations  were  used.   In  region 

IV  a  polynomial  curve  fit  was  used.   The  equations  defined  for 

each  region  are  listed  below. 

Region  I  -  (  0-18  knots) 

FSD  =  CDX1  (  u/  umax  )3 
Region  II  -  (18-22  knots) 

FSD  =  CONSTANT   -  CDX-  u 
Region  III  -  (22-36  knots) 

FSD  =  CDX2  u  -  CONSTANT 
Region  IV  -  (above  36  knots) 

FSD  =  CM3  u2  +  C^/u1'5 

3 .   Identification  of  Surge  Drag  Coefficients 

In  region  I,  C   _  is  equal  to  the  first  peak  value 

of  drag.   In  regions  II  and  III,  C^,.,.  and  CDX2  are  equal  to 

the  slopes  of  the  drag  curve  and  the  constants  were  solved 

for  by  applying  the  boundary  conditions  at  the  intersections 

of  the  adjacent  regions.   The  values  of  C^-  and  Cpv*  were 

computed  by  taking  F.SD  at  two  different  known  values  of  u 

and  solving  for  C„v,  and  C~v,  simultaneously. 
to      DX3      DX4 

B.   LINEARIZED  MODEL  ANALYSIS 

Since  the  surge  model  of  the  RTS5D  is  extremely  non- 
linear over  the  complete  speed  range,  analysis  is  very  com- 
plex.  Therefore  the  model  was  linearized  and  it's  performance 
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analyzed  at  various  points  on  the  FSD  curve.   Linearization 
of  the  model  was  accomplished  by  computing  the  linear 
coefficients  through  partial  differention  with  respect  to 
surge  in  each  of  the  four  region  equations  and  evaluating 
these  terms  at  various  operating  point  surge  velocities  (u  ) 
The  linearized  model  is  shown  in  Fig.  6  and  the  linear  coef- 
ficients for  the  surge  drag  force  FSD  are  shown  below. 


Region 

I     FSD,..  x  =  1  C™  (   u  /  u    ^2 


(■>,    >   =  3  GnY1  (  u  /  u   )' 
luQ;      DX1    o'   maxy 


ii   fsd'(Uo)  =  -  cDX5 

III    FSd'(   ,  =  CDX2 


IV    FSD(   .  =  2  CDX3  uQ  -  1.5  CDxi/  uQ2-5 


A  map  of  the  linearized  surge  directional  roots  are 
shown  in  Fig.  7.  The  characteristic  roots  for  each  operating 
region  are  shown  coded  [see  legend) .  Note  that  only  Region  I 
and  Region  IV  roots  are  u  dependent  and  therefore  move  as  uq 
changes.  The  direction  of  change  along  the  real  a_xis  is 
noted  by  the  arrows  shown.  Region  II  and  III  roots  are  fixed 
in  locations  and  independent  of  u  . 
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Figure  6 
Linearized  Surge  Model 

Region  IV  of  Fig.  7  is  of  particular  interest  since  the 
ability  to  operate  in  this  speed  range  is  a  primary  consi- 
deration in  the  development  of  the  SES.   It  can  be  seen 
from  the  map  of  the  roots  that  within  this  region,  at  speeds 
from  36  to  almost  60  knots,  poles  exist  in  the  RHP.   This 
is  a  result  of  the  negative  drag  slope  where  the  retarding 
force  actually  decreases  with  an  increase  in  speed,  the 
magnitude  of  variation  being  dependent  on  sea  state.   Region 
II  also  exhibits  this  negative  drag  characteristic  but  since 
it  is  of  such  a  narrow  speed  range.   operation  in  this  re- 
gion was  not  considered  in  this  design.   Further  modeling 
changes  for  the  other  degrees  of  freedom  were  based  on 
Region  IV  response.   In  Regions  I,  III  and  in  Region  IV  above 
60  knots,  the  system  will  be  stable  at  the  speed  u  where  FSD 
is  equal  to  commanded  thrust. 
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Figure  7 
Linearized  Surge  Root  Map 
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If  it  were  necessary  to  operate  at  a  speed  within  the 
negative  drag  slope  characteristic,  it  would  be  necessary  to 
actively  modulate  the  thrust  command  in  order  to  maintain 
a  nearly  constant  speed,  but  this  would  require  a  great  deal 
of  operator  attention  since  the  negative  drag  slope  in  Region 
IV  is  one-  of  slowly  divergent  instability.   For  this  reason 
a  closed  loop  speed  controller  has  been  included  in  the  data 
base  model  and  therefore  was  also  added  to  the  RTS-5D  model. 

C.   SIMPLIFIED  PROPULSION:  DYNAMICS 

Before  adding  a  speed  autopilot  to  the  model  one  addi- 
tional refinement  was  included  in  the  design.   Thrust  command 
changes  in  the  RTSSD  as  reported  in  Ref.  I  were  instantan- 
eously entered  into  the  equations  of  motion  as  thrust 
changes.   This  was  not  the  case  for  the  data  base  simulation 
program.   Specific  thrust  levels  are  generated  by  specific 
water  flow  rates  through  water  jet  nozzles.   These  flow 
rates  are  nearly  proportional  to  the  pump  speed  that  pro- 
duces them.   Since  the  pumps  are  driven  directly  by  gas 
turbine  engines,  a  change  in  t-hrust  occurs  with  a  change  in 
turbine  speed.   A  thrust  command  translates  to  a  power  tur- 
bine speed  command  which  as  an  output  has  power  turbine 
speed  achieved.   This  output  is  translated  into  thrust 
achieved.   It  was  specified  in  Ref.  7  that  one  of  the  goals 
in  the  base  line  design  was  to  have  power  turbine  speed 
control  response  characteristics  which  are  relatively  inde- 
pendent of  the  power  setting.   Accordingly,  in  the  design 
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of  the  auto  pilot  of  Ref.  7,  the  simplified  propulsion 
block  shown  in  Fig.  8  was  used.   The  magnitude  of  the  pro- 
pulsion system  time  constant,  x  was  not  specified.   For 
the  purpose  of  this  analysis  x  was  chosen  to  be  1  second. 


S  Fuel 
Command 


l  +Ts 


JT 


Figure  8 
Simplified  Propulsion  System 


D.   ANALYSIS  OF  LINEARIZED  MODEL  WITH  SPEED  CONTROLLER 

The  complete  linearized  model  with  speed  controller  that 
is  included  in  the  model  for  the  purposes  of  this  study  is 
shown  in  Fig.  9.   The  saturation  effects  of  fuel  control  are 
ign.ored  and  the  linear  representation  is  shown  as  K  which 
is  assumed  to  be  unity.   The  controller  includes  propor- 
tional plus  integral  control.   The  proportional  control  was 
included  to  keep  the  ship  speed  error  low  while  t.he  integral 
control  was  included  in  order  to  z:ero  out  long  term  control 
error . 
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The  characteristic  equation  for  the  system  in  Fig.  9  is 

S3  +  (1  +  FSD(U  \/   m  )  S2  +  l/m  (FSd',  )   +  K)  S  +  K/  mf  =0 

o  *  o 

Using  the  Routh  Hurwitz  Criterion  for  stability 

s3  1  l/m  (FSD'(u  j  +  K) 

S2         <  1  +  FSD'/   y/  m  )  K/mf 

o; 

S1   l/m  (1  +  PSD/   >/m)  (FSd',   %  +  K)  -  K/mT        0 

v  o;        K    o' 

s0  K/  m  r  0 

For  Stability  K  must  be  greater  than  0  and  greater  than 

FSDu0)  r  +  mrFSD'(uo) 
FSD(u  )r  +  »r.   m 

o 

The  values  of  controller  parameters  K  and  x  can  be  initialized 
by  the  operator  for  either  loosely  controlled  speed  or  tightly 
controlled  speed  for  such  situations  as  station  keeping, 
UNREPS,  or  test  operations.   For  comparison  to  the  data  base, 
these  parameters  were  selected  to  achieve  closest  agreement 
between  the  values  of  surge  achieved  by  the  RTS5D  and  DBSIM5D 
models  after  the  completion  of  a  360°  turn  with  50  knot  initial 
surge  veloctiy  and  a  thrust  effector  angle  of  15°.   This  test 
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condition  for  the  nonlinear  models  was  chosen  for  two  reasons 
First,  given  fuel  conservation  as  a  necessity,  50  knots  is  a 
likely  high  speed  initial  condition  because  it  is  close  to 
minimum  drag.   Second,  this  initial  condition  provided  an 
opportunity  to  test  the  operation  of  the  speed  controller  in 
the  negative  drag  region  of  FSD  since  the  decrease  in  speed 
going  into  the  turn  with  u  =50  knots  caused  increased  drag. 
The  values  of  controller  parameters  that  caused  closest  agree 
ment  between  the  two  nonlinear  models,  when  tested  in  the 
linearized  model  showed  instability  in  the  linear  sense. 
This  is  desirable  in  a  turning  mode  for  the  nonlinear  model 
because  it  allows  speed  to  decrease  and  the  turn  to  be  com- 
pleted more  quickly.   Stiff  speed  control  could  be  achieved 
by  increasing  the  controller  gains,  but  would  normally  be 
used  only  in  straight  ahead  runs. 

Coupling  of  the  equations  of  motion  developed  in  Chapter 
II  necessitated  a  repetition  of  the  parameter  identification 
procedure  resulting  in  new  values  for  C„y,  £w,  C_p  and  C„„p. 
These  values  were  adjusted  for  the  best  curve  fit  compared 
to  the  data  base  using  the  linearized  model  for  the  RTS5D 
in  a  360°  turn  using  15°  thrust  effector  angle  with  initial 
surge  velocity  50  knots.   An  identical  series  of  tests  to 
those  in  Chapter  III  were  conducted  comparing  the  new  model 
of  the  RTS5D  to  the  DBSIM5D.   Complete  validation  results 
are  tabulated  in  Chapter  VI. 
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E.  EFFECTIVE  THRUST  EFFECTOR  ANGLE 

As  a  result  of  the  negative  drag  characteristics  imposed 
on  the  model ,  RTS5D  yaw  rate  (r)  was  high  in  the  low  speed 
test  and  low  in  the  high  speed  test.   This  opposite  error 
at  both  ends  of  the  test  speed  spectrum  was  reduced  by 
introducing  a  correction  factor  to  the  thrust  effector  angle 
proportional  to  yaw  rate 

Z  --  *  Z  -  SLIP  •  r 
ef  f 

and  adjusting  2w  to  increase  r.   This  had  the  desired  effect 
of  reducing  error  at  both  ends  of  the  test  spectrum  since 
the  reduction  of  the  thrust  effector  angle  was  greater  at 
the  high  end. 

F.  BROACH  CONDITION  FLAG 

As  stated  in  Chapter  I,  broach  response  characteristics 
are  not  included  in  this  model,  therefore  a  broach  condition 
warning  flag  was  included  in  the  design  in  order  to  alert 
the  operator  who  could  then  take  corrective  action  to  avoid 
broaching. 

Broaching  is  caused  by  a  water  intake  being  exposed  to 
air  and  is  a  function  of  surge  speed,  roll  angle,  pitch 
angle,  drift  angle  and  sea  state.   For  the  data  base  3K-SES 
with  broach  boundaries  defined  at  10%  air  by  volume,  the 
broaching  boundaries  in  calm  water  in  terms  of  roll  angle, 
pitch  angle  and  drift  angle  are  shown  in  Fig.  11  for  40  knots 
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and  Fig.  12  for  60  knots.   Reference  6  listed  maximum  thrust 
effector  angles  as  a  function  of  surge  speed  and  air  flow 
in  order  to  avoid  broaching.   This  condition  is  shown  plotted 
in  Fig.  10  for  fixed  plenum  air  flow  rate. 


EFFECTOR 

ANGLE 

(DEG) 


BROACH 


0   (KNOTS) 


SURGE 
VELOCITY 

Figure  10 
Maximum  Thrust  Effector  Angle  As  A  Function 
Of  Surge  Velocity  In  Order  To  Avoid  Broaching 

As  stated  in  Chapter  III  the  model  developed  in  Ref.  1  did 
not  flag  the  broaching  condition.   To  implement  the  broach 
warning  flag  and  display  the  other  new  parameters  for  the 
modified  RTS5D  the  pilot  graphic  display  described  in  Ref.  1 
was  changed  to  incorporate  the  necessary  data.   Figure  13A 
depicts  typical  information  displayed  to  the  pilot  under 
normal  operating  conditions.   The  new  information  to  the 
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operator  shown  in  the  lower  right  of  Fig.  13A  is  sea  state 
(SEAST) ,  computer  assist  (CASST,  on=l  off=0) ,  controller 
gains  (K,KK)  and  speed  commanded  (SPCOM) .   If  the  operator 
exceeds  the  thrust  effector  limits  (ZMAX)  to  prevent 
broaching  (see  Fig.  10),  a  flashing  warning  is  displayed 
in  place  of  the  bottom  line  of  information  on  the  screen. 
This  is  shown  in  Fig.  13B. 
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VARIABLE  DYNAMICS 

U    (HEAVY)  r^ 
UDOT( LIGHT) 


SEA  TRACK 


V    (HEAVY) 
VDOT( LIGHT) 

R    (HEAVY) 
RDOT( LIGHT) 

P    (HEAVY) 
PDOT (LIGHT) 

Q    (HEAVY) 

QDOT( LIGHT) 


S^ 


nr* 
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-J=-mA 


BW* 
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CONTROL  INPUTS     TIME  NAVIGATIONAL  DATA 

THRUST   RUDDER  DELT  ELAPSED  PITCH  HEADING/RATE  DRIFT  ROLL  SPEED 

300000   -.20   .084    22.2  .52    325/2.1    -.28    .8   50.1 

T7      T8      T9      T10  SEAST  CASST   K    KK    SPCOM 


75000   75000   75000 


750000 


2900  3 
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Figure  13A 
Pilot  Graphic  Display  (Normal) 
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VARIABLE  DYNAMICS 


U        (HEAVY) 
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V        (HEAVY) 
VDOT( LIGHT) 

R        (HEAVY) 
RDOT(LIGHT) 

P        (HEAVY) 
PDOT( LIGHT) 

Q        (HEAVY) 

QDOT(LIGHT) 
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SEA   TRACK 
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CONTROL  INPUTS     TIME  NAVIGATIONAL  DATA 

THRUST   RUDDER  DELT  ELAPSED   PITCH  HEADING/RATE  DRIFT  ROLL  SPEED 

3000000   -.52    .084    22.2    .$2  325/2.1    -.28    .8   50 . 1 

T7      T8      T9      T10     SEAST  CASST   K    KK    SPCOM 

—  B   R   0   A   C   H  — 

Figure  13B 
Pilot  Graphic  Display  (Broach  Condition  Exists) 
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V.   IMPLEMENTATION 

A.  INTRODUCTION 

Given  the  existence  of  a  Real  Time  Simulation  of  a  SES , 
a  decision  was  made  to  refine  that  model  using  more  accurate 
equations  of  motion,  improved  computer  iteration  time  and 
operator  hardware  modifications.   A  complete  description  of 
the  RTS5D  is  contained  in  Ref.  1,  therefore  only  changes 
are  discussed  here,  and  included  are  the  introduction  of 
negative  drag  in  a  sea  state,  propulsion  dynamics,  speed 
control  and  a  broach  condition  warning  in  software.   Hard- 
ware modification  included  changing  the  thruster  control 
box  to  make  it  more  accessable  to  an  operator. 

B.  REQUIREMENTS 

The  following  criteria  were  required  of  the  remodeled 
RTS5D. 

1)  The  refined  simplified  equations  of  motion  would  be 
solved  and  the  results  output  on  a  real  time  basis. 

2)  The  solution  would  be  subject  to  real  time  control 
efforts  generated  by  an  operator  observing  the  output. 

3)  Computer  iteration  time  was  required  to  be  less  t-han 
laQi  ms. 

A  more  complex  discussion  of  condition  3)  is  presented  in 
Appendix  A. 
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C.  HARDWARE  DESCRIPTION 

The  hardware  arrangement,  in  block  diagram  form  is  shown 
in  Fig.  14.   The  only  change  to  the  RTS5 D  was  the  redesign 
and  relocation  of  the  thruster  console  to  the  pilots  seat 
where  it  can  be  more  conveniently  utilized  by  the  operator. 
The  thruster  console  uses  "linear  movement"  potentiometers 
in  place  of  previously  used  "dial"  potentiometers  in  order 
to  more  realistically  simulate  thrust  controls  for  the  SES. 
The  thruster  console  is  shown  in:  Fig.  15. 

D.  SOFTWARE  DESCRIPTION 

The  remodeled  software  package  was  a  FORTRAN  IV  digital 
program  which  consisted  of  a  main  graphics  program  and  three 
major  subroutines.   It  was  necessary  to  make  a  third  sub- 
routine out  of  the  equations  of  motion  computation  that  had 
previously  resided  in  the  main  program  prior  to  making  any 
changes  to  these  equations.   It  was  discovered  when  the 
first  program  modifications  were  attempted  that  the  compiler 
for  the  XDS-9300  computer  was  operating  at  its  limit  with 
the  program  as  it  then  existed,  therefore,  before  changes  to 
the  equations  could  be  implemented,  the  program  structure 
had  to  be  modified.   The  failure  mode  was  eliminated  from 
the  new  model  because  of  initial  assumptions  that  the  thrusts 
would  be  all  applied  at  the  same  point.   Thrust  failures 
should  not  be  considered  until  the  model  includes  broaching 
dynamics  which  is  a  major  source  of  thrust  loss.   The  flow 
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chart  for  the  remodeled  RTSBD  is  shown  in  Fig.  Id.   The  com- 
plete multiplexing  algorithm  is  shown  in  a  flow  chart  in 
Fig.  17. 


56 


XDS 
9300 

DIGITAL 

COMPUTE!* 


COMCOR 
CI  5000  ANALOG 
COMPUTER 


-FT. 

PLOTTER 
8  CH 
STRIP 
CHART 
REC 


FT 


TTY-35 
TELETYPE- 
WRITER 


LINE 
PRINTER 


DISC 
DRIVE 


1 


DISC 
DRIVE 


/ANALYST 
I  SEAT 


THRUSTER 
CONSOLE 


b 


JOYSTICK 


ADAGE 
AGT-10 


GRAPHICS 

DISPLAY 

CONSOLE 


GRAPHICS! 

DISPLAY 

CONSOLE 


ADAGE 
AGT-10 
5F 


Figure  Ik 
RTS5D  Mod  Block  Diagram 
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Figure  15 
Thruster  Console 
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Figure  16 
RTS5D  Mod  Program  Flow  Chart 
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Figure  17 
RTS5D  MOD  DISPLAY  MULTIPLEX  ALGORITHM 
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VI.   RTS5D  MODS  I  AND  II  RESPONSE  CHARACTERISTICS 

RTS5D  MOD  I  includes  negative  drag  effects  in  a  sea  state, 
simplified  propulsion  dynamics,  speed  control  and  a  broach 
condition  warning  flag.   RTS5D  MOD  II  is  identical  to  MOD  I 
except  that  it  includes  the  correction  factor  for  thrust 
effector  angle  developed  in  Chapter  IV,  Section  E. 

The  validation  of  the  modified  RTS5D  was  accomplished 
in  the  same  manner  as  the  RTS5D.   CSee  Chapter  III)  ..   The 
results  are  displayed  in  the  following  tables.   A  complete 
360°  turn  comparison  of  the  DBSIM5D,  the  RTS5  D  and  the  modi- 
fied RTS5 D  at  50  knots,  15°  thruster  angle  in  shown  in  Fig. 
18.   Additionally,  a  comparison  of  transient  response  was 
undertaken  and  the  time  of  first  peak  for  the  variable  v,  r, 
4>  and  9  for  each  model  was  graphed  in  Figs.  19-2  1  for  speeds 
of  40.  ,  5a  and  60  knots.   Transient  response  characteristics 
are  an  important  consideration  in  the  development  of  a  real 
time  model  and  are  discussed  in  Appendix  A.   The  transient 
response  of  the  RTS5~D  models  was  faster  than  that  of  the 
data  base  model  with  the  exception  of  roll  angle  O) .   This 
impacted  on  the  validity  of  a  real  time  solution  of  the  equa- 
tions of  motion  and  the  display  of  the  results. 
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TABLE  IV 

RTS5D  Mod  I 

and 

Mod 

II 

Performance 

Test  at 

^0 

Knot 

s 

effector 
angle 

RTS5D 

Mod  : 

[ 

RTS5D  Mod 

II 

4> 

in  degrees 

<J> 

5° 

1st  pk 

.52 

Qss 
.45 

1st  pk 
.40 

Qss 
.38 

10° 

.96 

9 

.92 
in  degrees 

.79 

6 

.76 

5° 

1st  pk 
1.21 

Qss 
1.17 

1st  pk 
1.20 

Qss 
1.18 

10° 

1.20 

u 

1.18 
in  ft/sec 

1.20 
u 

1.18 

5° 

1st  pk 
62.56 

Qss 
62.64 

1st  pk 
62.58 

Qss 
67.65 

10° 

n/a 

V 

61.37 
in  ft/sec 

n/a 

V 

61.71 

5° 

1st  pk 
2.62 

Qss 

2.47 

1st  pk 
2.33 

Qss 

2.27 

10° 

3.53 

r 

3.52 
in  degrees 

3.23 
/sec   r 

3.20 

5° 

1st  pk 
.69 

Qss 
.59 

1st  pk 

.57 

Qss 
.54 

10° 

1.18 

1.22 

1.04 

1.09 
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TABLE  V 
RTS5D  Mod  I  and  Mod  II  Performance  Test  at  50  Knots 

RTS5D  Mod  I  RTS5D  Mod  II 


effector 

angle 

5° 

1st  pk 
.50 

10° 

.90 

15° 

1.31 

5° 

1st  pk 
1.23 

10° 

1.24 

15° 

1.24 

5° 

1st  pk 
n/a 

10° 

n/a 

15° 

n/a 

5° 

1st  pk 

2.  57 

10° 

3.48 

15° 

4.13 

5° 

1st  pk 

.57 

10° 

.96 

15° 

1.31 

u 


in 

degrees 

* 

Qss 
.39 

1st  pk 

.43 

Qss 
.39 

.82 

.83 

.78 

1.36 

1.24 

1.19 

in 

deg: 

rees 

9 

Qss 
1.19 

1st  pk 
1.23 

Qss 
1.20 

1.19 

1.23 

1.20 

1.18 

1.23 

1.19 

in 

ft/ 

sec 

u 

Qss 

82.73 

1st  pk 
n/a 

Qss 
82.98 

78.68 

n/a 

79.55 

68.30 

n/a 

72.75 

in 

ft/ 

sec 

V 

Qss 
2.30 

1st  pk 

2.39 

Qss 
2.29 

3.33 

3.25 

3.25 

4.29 

3.98 

4.00 

in 

deg 

rees/sec 

r 

Qss 
.39 

1st  pk 
.49 

Qss 
.41 

.85 

.88 

.86 

1.62 

1.25 

1.41 
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TABLE  VI 

RTS5D  Mod  I  and  Mod  II  Performance  Test  at  6  0  Knots 

effector 
angle  RTS5D  Mod  1  RTS5D  Mod  II 


ko 


<}>         in  degrees  $ 

1st  pk  Qss  1st  pk        Qss 

.53  .39  .46  .41 

9         in  degrees  9 

1st  pk  Qss  1st  pk        Qss 

1.23  1.20         1.24  1.20 

u  u 

1st  pk  Qss  1st  pk        Qss 

n/a  100.21       n/a  99.59 

v  v 

1st  pk  Qss  1st  pk        Qss 

2.64  2.29         2.49  2.35 


1st  pk  Qss  1st  pk        Qss 

5°         .52  .32  .46  .36 
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360  Turn  Comparison  of  Four  Models 
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Response  Time  Comparison  at  ^0  Knots 
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Figure  21 
Response  Time  Comparison  at  60  Knots 
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VII.   CONCLUSIONS 

The  RTS5D  has  been  shown  in  Section  VI  to  be  a  viable 
input/output  model  of  a  3K-SES  by  closely  matching  its  output 
characteristics  for  a  given  input  to  those  of  the  data  base 
model. 

The  refinement  of  the  right  side  of  the  simplified  equa- 
tions of  motion,  especially  the  inclusion  of  the  negative 
drag  characteristic  in  a  sea  state  has  made  the  model  more 
accurate  over  a  wide  speed  range  and  gives  the  operator,  upon 
program  initialization,  the  ability  to  choose  different  drag 
characteristics  caused  by  sea  state. 

The  model  was  tested  in  a  negative  drag  region  with 
automatic  speed  control  and  simplified  propulsion  dynamics 
and  it  was  demonstrated  that  by  properly  adjusting  controller 
gains  that  the  output  characteristics  of  the  DBS IMS Q  were 
closely  matched. 

Previously  conducted  thrust  failure  analysis  with  man 
generated  corrective  action  has  been  shown  to  be  invalid 
until  broach  dynamics  are  included  into  the  model.   Previous 
tests  allowed  maximum  thrust  effector  angle  correction  regard- 
less of  this  consideration.   The  modified  version  of  the  RTS5  D 
does  not  include  complete  broach  dynamics  but  does  warn  the 
operator,  if  he  should  exceed  a  maximum  thrust  effector  angle 
for  a  given  forward  speed. 
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The  3K-SES  may  be  operated  in  the  negative  drag  slope 
speed  region  approximately  C 36 -5  8)  knots  without  the  aid  of 
a  speed  controller  but  would  require  constant  manual  modula- 
tion of  the  thrusters.   Since  minimum  manning  is  one  of  the 
criteria  in  3K-SES  design:,  a  speed  controller  is  necessary 
to  free  operator  attention  to  other  duties  when  operating 
in  this  speed  range. 

AnotFter  consideration  of  the  implementation  of  the  non- 
linear drag  curves  is  that  for  a  given  sea  state,  two  values 
of  minimum  drag  can  be  determined,  one  in  Region  IV  and  one 
"between  Region  II  and  Region  III.  Knowledge  that  drag  can 
be  lower  in  a  particular  operating  region  for  fuel  conserva- 
tion is  of  prime  importance. 

It  is  recognized  that  a  time  delay  exists  in  the  solution 
of  the  equations  of  motion  and  a  time  delay  exists  in  the 
graphical  display  of  the  solutions  to  the  operator.   Real 
time  interaction  between  the  man  in  the  loop  and  the  computer 
is  achieved  by  setting  the  integration  step  size  equal  to 
th_e  computer  loop  iteration  time.   Real,  time  consi  deration-S 
are  further  discussed  in  Appendix  A. 


70 


VIII.   RECOMMENDATIONS 

The  following  recommendations  are  made  concerning  future 
improvements  in  the  RTS5D  model. 

1.  Further  refinement  of  the  right  hand  side  of  the 
equations  of  motion  is  necessary  in  order  to  more  closely 
match  the  characteristics  of  the  data  base  model. 

2 .  Implement  software  broach  dynamics  in  the  model  with 
thrust  applied  differentially  at  different  lever  arms  Csee 
Fig.  3 )  so  that  thrust  failures  and  effectiveness  of  correc- 
tive action  may  be  evaluated. 

3.  Rewrite  the  program  such  that  the  AGT-10's  would 
provide  asyncronous  parallel  graphics  C existing  equipment) 
or  implement  the  model  on  a  more  modern  system  which  uses 
parallel  processing  and  improved  graphics. 

■+.   Introduce  the  sixth  degree  of  freedom  Oheave}  as  a 
parallel  operation.   Since  heave  time  constant  is  much  faster 
than  the  rest  of  the  system  dynamics,  real  time  heave  analysis 
could  be  conducted  simultaneously  if  cross  coupling  is  neg- 
lected. 

5  .   The  system  is  designed  for  man  in  the  loop  interface, 
so  further  operator  te-sting  is  recommended  as  a  valuable 
source  of  information  for  future  refinement  of  the  model. 
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APPENDIX  A 
REAL  TIME  ANALYSIS 

T.  S.  Nelson,  in  Ref.  1,  accomplished  real  time  solution 
of  SES  equations  of  motion  with  the  use  of  a  1000  Hz  clock 
to  measure  computer  iteration  time  and  by  using  that  time  as 
an  integration  step  size  for  rectangular  integration  to  solve 
the  simplified  dynamic  equations.   He  concluded,  by  comparing 
the  solutions  of  these  equations  to  that  obtained  by  using  a 
Runga  Kutta  variable  integration  step  size,  that  the  computer 
iteration  loop  time  must  stay  below  100  ms  in  order  to  main- 
tain at  least  3%  accuracy  of  solution.   He  further  developed 
a  multiplexing  algorithm  in  order  to  output  these  solutions 
to  an  operator  in  the  form  of  graphics  on  a  real  time  basis. 
This  model  can  be  invalidated  as  a  real  time  model  in  two 
ways.   First,  in  order  to  have  real  time,  the  integration 
step  sizie  used  must  be  equal  to  the  computer  iteration  time 
and  in  this  model,  because  of  the  multiplexing  algorithm, 
each  loop  time  is  different.   Th-e  integration  step  size, 
dt  determined  by  measuring  loop  N  is  used  to  solve  the  equa- 
tions in  loop  N+l,  which  has  a  different  loop  time,  there- 
fore the  solutions  are  not  exactly  real  time.   Secondly, 
the  multiplexing  scheme  makes  the  model  into  a  sampled  data 
system  with  sample  time  (Tsl  approximately  Z  seconds.   Even 
if  the  first  condition  did  not  apply  and  real  time  solutions 
to  the  equations  could  be  assumed,  then  the  output  to  the 
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operator  would  be  sampled  real  time  data.   Given  the  fastest 
time  constant  for  the  RTS5D  simplified  equations  (approxi- 
mately 1.3  seconds  for  pitch  rotational  motion),  it  can  be 
seen  that  the  Nyquist  Criteria  is  clearly  violated.   Yet 
even  with  these  obvious  limitations,  it  has  been  demonstrated 
that  the  RTS5D  can  closely  reproduce  the  output  characteris- 
tics of  the  DBSIM5D  and  therefore  is  a  useful  model.   Surge 
speed  as  well  as  sway,  yaw  and  roll  (Figs.  21-2  3)  have  much 
longer  effective  time  constants.   It  is  recognized  that 
there  is  a  delay  between  the  real  time  of  any  event  and  the 
observation  of  that  event,  therefore  "Real  Time"  as  used  in 
this  report  is  a  matter  of  degree,  or  more  specifically 
how  much  delay  time  between  solution  and  display  can  be 
tolerated.   Given  the  scale  of  graphics  presented  to  the 
operator  in  this  model  and  the  speed  with  which  the  SES 
responds,  the  RTS5 D  is  an  acceptable  model,  however  it 

would  be  desirable  to  reduce  the  sample  time  T   and  make  the 

s 

integration  step  size  equal  to  the  computer  loop  iteration 
time  for  each  iteration. 

An  algorithm  to  set  the  integration  step  size  equal  to 
the  computer  iteration  time  of  the  current  loop  was  imple- 
mented as  a  software  design,  oh.  ange  to  the  RTS5  D  model.   This 
was  accomplish-ed  by  defining  the  integration  step  size 
as  an  array  DELTC  N)  with  N  equal  to  the  total  number  of  itera- 
tion's ne-cessary -t  o  display  all  of  t  he  graphical  information 
to  the  operator.   Each  loop  iteration  time  was  ne  asured  and 
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stored  in.  this  array  during  the  first  pass  through,  the 
multiplexing  algorithm.   DELTCMJ  was  initialized  at  z.ero 
so  that  time  would  not  increment  and  integration  would  not 
be  performed  during  the  first  pass  through  the  multiplexing 
algorithm.   On  subsequent  passes  through  the  multiplexing 
algorithm  the  integration  step  size  was  equal  to  the  loop 
time  in  which  it  was  used. 

"This  timing  algorithm  when  tested  in  the  RTS5  D  model 
was  found  to  accurately  represent  real  time  but  the  devia- 
tion from  previously  used  integration  step  size  values 
caused  instability  in  both  the  pitch  and  roll  dynamics. 
This  needs  further  investigation. 

Sturgeon,  in  Rex.  8,  demonstrated  the  value-  of  using 
parallel  array  processors  to  reduce  computational  time  in 
real  time  simulations.   Such  a  model  could  be  designed  for 
the  RTSSD  using  the  existing  equipment.   The  AGT-10's  which 
are-  used  for  displaying  the  output  can  inte-rface  with  the 
XDS-93Q0  main  computer  in  two  ways.   The  multiplexing  scheme 
of  the  RTS5D  uses  a  "hand  shaking"  method  with  a  GRAPHO  or 
TEXTO  subroutine  call  in  each  loop.   It  is  precisely  the 
inclusion  of  graphics  in  each  loop  that  causes  problems  for 
real  time  presentation.   The  solution  of  the  equations  requires 
only  approximate- ly  3 5  ms .   Each  GRAPHO  or  TEXTO  call,  with 
its  associated  "hand  shaking"  time,  makes  up  the  rest  of  the 
loop  time,  approximately  4  0  ms  more.   It  can  be  seen  that  by 
introducing  two  such  calls  in  one  loop  that  the  10  0  ms  time 
restriction  would  be  exceeded. 
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An  improvement  to  this  method  would  be  to  use  the  second 
interface  between  the  AGT-10's  and  the  XDS-9300  which  is 
that  the  AGT-10's  have  a  direct  connection  to  some  memory 
locations  in  XDS-9300.   By  solving  the  dynamics  for  th,e  SES 
and  storing  the  data  in  these  memory  addresses  in  the  form 
of  arrays,  the  AGT-TTrs  could  be  made  to  operate  asyncronously 
in  parallel  as  array  processors  by  using  a  process  known  as 
"cycle  stealing"  to  fetch  the  data  from  the  XDS-9300  and  to 
display  this  on  the  screen  for  the  operator.   This  process 
would  have  the  advantage  of  reducing  and  making  more  uniform 
the  computer  iteration  time  used  as  the  step  size  for  the 
rectangular  integration  in  the  solution  of  the  SES  equations 
of  motion.   The  sample  time  would  also  be  reduced  because  all 
of  the  time  wasted  in  the  "hand  shaking"  is  eliminate-d.   It's 
drawback  is  that  the  output  is  not  truely  real  time  because 
of  the  asyncronous  operation  of  the  AGT-10's.   ffowever,  the 
reduced  sampled  data  time  and  computer  iteration  time  would 
make  it  a  closer  approximation  of  real  time  than  presently 
exists . 
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APPENDIX  E 
RTS5D  MODIFIED  PROGRAM  NOMENCLATURE 


All 
A21 
A22 

A31 

A33 
A34 

AM_ 

A3SYD 

A23Y0D 

AD 

A1GOA 

AM 

APITCH 

APRINT 

ARR 

A  RATE 

ARCLL 

ASTOP 

AXST 

AYST 

2ETA 

BOAT 


A/D  trunk  500  voltage 

A/D  trunk  501  voltage 

Added  mass  coefficient  in  yaw  moment  equation 

A/D  trunk  502  voltage 

Added  mass  coefficient  in  roll  moment  equation 

Added  mass  coefficient  in  pitch  moment  equation 

A/D  trunk  503  voltage 


array  of  A/D  lines 

speed  line  index  counter 

mass 

pitch  angle  in  degrees 

;f  iterations  between  print  execute  commands 

Real  IARR  (l) 

turn  rate  in  degrees 

roll  angle 

performance  index  J\  limit 

initial  condition  X-coordinate  on  sea  track 

initial  condition  Y-coordinate  on  sea  track 

drift  angle 

subroutine  to  generate  perspective  of  SSS 
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GA33T 

CDP 

CDX(., 

GDY 

CDZP 

DAL 

DELT 

DELTA 

DFCOM 

DFCOM  1 

DFCOM  2 

DTI  ME 

DTIMPLT 

DUG  CM 

DXPLCT 

DYPLCT 

DYREPET 

EFLG 

ETIME 

EYRE PET 

EXPLOT 

EYPLCT 

Fl 

F2 

F3 


automatic  speed  control  mode  indicator 

bow  real  forces  lumped  coefficient 

drag  forces  lumped  coefficients 

sway  forces  lumped  coefficient 

sidewall  rolling  moment  lumped  ceofficient 

digital  to  analog  call 

iteration  loop  time 

commanded  iteration  loop  time 

total  change  in  fuel  command  from  speed  comtroller 

change  in  fuel  command  due  to  proportional  speed 

controller 

change  in  fuel  command  due  to  proportional  speed 

controller 

v  time 

scaled  value  for  v 

difference  between  ucom  and  u 

same  as  DTIMPLT 

scaled  value  for  v 

restart  value  for  v 

error  flag  for  TEXTC  call 

r  time 

restart  value  for  r 

scaled  r  time 

scaled  r 

stern  buoyancy  force 

bow  seal  pitch  force 

bow  buoyance  force 
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FCOM  fuel  command 

FSP  port  sidewall  buoyancy  force 

FSS  starboard  sidewall  buoyancy  force 

FTIME  v  time 

FTIMPLT  scaled  v  time 

FXPLOT  same  as  FTir.TLT 

FYPLCT  scaled  v 

FYREPET  restart  value  for  v 

G  gravitational  acceleration 

GRAPKO  subroutine  to  project  image  (non  alpha-numeric) 

GTIME  R  time 

GTIMPLT  scaled  R  time 

GXFLOT  same  asGTIMPLT 

GYPLCT  scaled  R 

GYREPET  restart  R  value 

H/.  >.  dynamics  maximum  and  minimum  array 

HEAD  craft  heading 

HOLD  subroutine  to  freeze  display/program 

HTIME  u  time 

HTIMPLT  scaled  u  time 

HYREPET  analog  restart  y  coordinate  on  sea  track 

HYSESX  analog  y  position  of  craft 

HYSESY  analog  x  position  of  craft 

IARR(l)  digitized  mode  selector 

IARR(2)  digitized  effector  angle 

IARR(3)  digitized  T? 
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IARR(4)  digitized  T8 

IARR(5)  digitized  T9 

IARR(6)  digitized  T10 

IARR(7)  digitized  surge  velocity  command 

ICO  multiplex  graphics  index 

ICCA  speed  line  index 

ICRAFT  vehicle  designator 

IER  error  flag 

IIPR  time  history  limit 

IIPRD  time  history  limit 

IIPRE  time  history  limit 

IIPRF  time  history  limit 

IIPRG  time  history  limit 

IIPRH  time  history  limit 

IIPRC  time  history  limit 

IIPRP  time  history  limit 

IIPRQ  time  history  limit 

IIPRR  time  history  limit 

IIPRS  time  history  limit 

IIX  print  counter 

IJ  craft  perspective  counter 

IJA  craft  perspective  counter 

IJAD  time  history  counter 

IJAE  time  history  counter 

IJAF  time  history  counter 

I JAG  time  history  counter 

IJAH  time  history  counter 
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IJAO  time  history  counter 

IJAP  time  history  counter 

IJAQ  time  history  counter 

IJAR  time  history  counter 

IJAS  time  history  counter 

IJD  time  history  counter 

IJE  time  history  counter 

IJF  time  history  counter 

IJG  time  history  counter 

IJH  time  history  counter 

IJO  time  history  counter 

UP  time  history  counter 

IJQ  time  history  counter 

IJR  time  history  counter 

IJS  time  history  counter 

IL  graphics  console  display  array 

I LAI  graphics  console  display  array 

ILA2  graphics  console  display  array 

ILA3  graphics  console  display  array 

ILA4  graphics  console  display  array 

ILA5  graphics  console  display  array 

ILA6  graphics  console  display  array 

ILA7  graphics  console  display  array 

ILA8  graphics  console  display  array 

ILA9  graphics  console  display  array 

ILA10  graphics  console  display  array 
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ILA11  graphics  console  display  array 

I LAI 2  graphics  console  display  array 

I LAI 3  graphics  console  display  array 

I LAI 4  graphics  console  display  array 

I LAIS  graphics  console  display  array 

I  LAIS"  graphics  console  display  array 

ILA17  graphics  console  display  array 

I LAI 8  graphics  console  display  array 

ILA19  graphics  console  display  array 

ILA2Q'  graphics  console  display  array 

ILA21  graphics  console  display  array 

ILAZ2  graphics  console  display  array 

ILA23  graphics  console  display  array 

ILA2<+  graphics  console  dispLay  array 

ILA2S  graphics  console  display  array 

ILA25  graphics  console  display  array 

ILA2  7  graphics  console  display  array 

ILA2  8  graphics  console  display  array 

ILA29  graph.ics  console  display  array 

ILAP  graphics  console  display  array 

IPACK  subroutine  to  load  display  arrays 

IPLOT  graphics  console  display  array 

IPR  time  history  counter  limit 

IPRD  time  history  counter  limit 


IPRE  time  history  counter  limit 

IPRF  time  history  counter  limit 

IPRG  time  history  counter  limit 

IPRH  time  history  counter  limit 

IPRO  time  history  counter  limit 

IPRP  time  history  counter  limit 

IPRQ  time  history  counter  limit 

IPRR  time  history  counter  limit 

IPRS  time  history  counter  limit 

ISWA  display  array 

ITEXT  subroutine  to  display  alpha-numeric  data 

IVEIWA  craft  perspective/north  road  array 

IVIEWAA  sea  track  array 

IVIEWB  sea  track  axis  array 

IVIEWC  u  time  history  array 

IVIEWD  v  time  history  array 

IVIEWE  r  time  history  array 

IVIEWF  v  time  history  array 

IVIEWG  r  time  history  array 

IVIEWH  u  time  history  array 

IVIEWO  p  time  history  array 

IVIEWP  q  time  history  array 

IVIEWQ  p  time  history  array 

IVIEWR  q  time  history  array 

IVIEWS  analog  sea  track  array 

IVIEWZ  time  history  axis  array 
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IX  x  axis  moment  of  inertia 

J  time  history  counter 

JD  time  history  counter 

JE  time  history  counter 

JF  time  histroy  counter 

JG  time  history  counter 

JH  time  history  counter 

JO  time  history  counter 

JP  time  history  counter 

JQ  time  history  counter 

JH  time  history  counter 

JS  time  history  counter 

IC  proportional  speed!  controller  gain 

KK  intregal  speed  controller  gain 

L3  lever  arm  for  buoyancy  force 

LD  draft  of  sidewall 

LP  side  thrust  lever  arm 

LZ  lever  arm  for  sway  drag  force 

N  clock  interrupt  count 

NAD  AD  truni  line  array 

ONEMOT  subroutine  to  calculate  SES  dynamics 

CO  lever  arm  for  thrust  side  component 

OTIME  p  time 

OTIMPLT  scaled  p  time 

CXPLCT  same  as  OTIMPLT 

OYPLCT  scaled  p 

CYREPST  restart  p  value 
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p 

roll  rate 

PB 

plenum  pressure 

PDCT 

roll  acceleration 

PDOTM 

max  roll  acceleration 

PHI 

roll  angle 

PMAX 

max  roll  rate 

PSI 

craft  heading 

PTI'tfE 

• 

q  time 

PTIMPLT 

• 

scaled  q  time 

PXPICT 

same  as  PTIMPLT 

FYPLCT 

• 

scaled  q 

PYRSPST 

■ 

restart  q  value 

Q 

pitch  rate 

QDCT 

pitch  acceleration 

QDCTM 

• 

max  q 

QIvIAX 

max  q 

QTir.E 

p  time 

QTIT;1FLT 

scaled  p  time 

QXPLCT 

same  as  QTIMPLT 

QYPLOT 

scaled  p  value 

QYREPET 

restart  p  value 

R 

craft  turn  rate  r 

RDCT 

» 

r 

RDCTI-T 

max  r 

R2ADCLOCK 

subroutine  to  sample 
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RESET 

RHO 

RMAX 

RTIME 

RTIMP1T 

RX 

RXPICT 

RY 

RYPLCT 

RYREPET 

RZ 

SSAST 

SESX 

SESY 

SFl 

SF14 

SLIP 

SPDLIN 

SPEED 

STARTGLCCK 

STIME 

STIMPLT 

SXPLCT 

SYPLCT 

SXREPET 

SYREPET 


subroutine  to  reset  analog  computer 

density  of  water 

maximum  yaw  rate  of  craft 

q  rate 

scaled  q  time 

h 

same  as  RTir.TPLT 

TY 

scaled  q  value 

restart  q  value 

Iz 

sea  state 

scaled  X/. 
o 

scaled  Y 
o 

scale  factor  in  time  history  arrays 

scale  factor  in  time  history  arrays 

Beta  influence  coefficient 

speedline  lower  screen  limit 

Craft's  velocity  in  knots 

subroutine  to  start  clock 

analog  time 

scaled  analog  time 

HYSESX 

HYSESY 

restart  value  for  HYSESX 

restart  value  for  HYSESY 
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T  total  thrust 

T10  thruster  number  k 

T10MAX  max  T10  value 

T7  thruster  number  1 

T7T-1AS  max  T7  value 

T8  thruster  number  2 

T8KAX  max  value  of  T8 

T9  thruster  number  3 

T9MAX  max  value  of  T9 

TEXTO  subroutine  to  display  alpha-numeric  data 

TDCT  T 

THETA  pitch  angle 

Til.!  u  time 

TIME  time 

TIMPLT  scaled  TIM 

TINT  real  clock  interrupt  frequency 

TITLE  program  name  display  array 

TITLED  program  name  display  array 

TITLS1  program  name  display  array 

TM  manual  thrust  command 

TriAX  max  T 

TSIDE  sum  of  side  thrust  components 

TYAW  sum  of  yaw  moment 

U  forward  velocity  (surge) 

UCOM  forward  velocity  (surge)  command 
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UDCT 


UDOTM 
UMAX 

u:iul 

V 

VCD 
VDOT 
VDOTM 

vmax 

vs 
w 

WDOT 

WE 

WRITECLGGK 

WW 

X 

xc 

XCDOT 
XDX 
XGX 
XPLOT 

XREPET 
XST 
YO 
YCDCT 


u 

max  u 

max  u 

scale  factor  in  speedline  generation 

side  velocity  (sway) 

thruster  console  pot  array 

• 

v 

max  v 

max  v 

total  velocity 

heave  velocity 

heave  acceleration 

width  of  bow  seal 

subroutine  to  assign  clock  interrupts  to 

variable 

lever  arm  for  sway  drag  forces 

Yo 

X 

o 

Y 

Ao 

boat  perspective  scale  factor 

boat  perspective  scale  factor 

display  array 

reset  value  for  SESX 

scale  factor  for  analog  plot 


y 
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YHIGH  point  D  vertical  position  cons  tan" 

YPLOT  display  array 

YREFET  reset  value  for  SESY 

-31  scale  factor  for  analog  plot 

Z  operational  effector  angle 

Z2  upper  effector  angle  limitation 

Z3  lower  effector  angle  limitation 

Z^  zero  angle  upper  limit 

Z5  zero   angle    lever   lir.it 

ZA3  scale  factor  for  Z10 

ZBRC  effector  "breach  iir.it  angle 

ZI  -edified  Z 

ZZZZ  fixed  step  effector  angle  input 
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Appendix  D 
RTS5D  Modified  Wiring  Diagram 
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